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The present invention relates to a optimising method for vibration damping treatments ol 
5 vehicle body panels according to the preamble of claim 1 . 

The current trend in the automotive industry towards reductions of treatment weight and 
cost, along with increasing demands for shorter development time and improved vehicle 
Noise and Vibration Harshness (NVH) characteristics, dictates the necessity for design 
10 methodologies which allow efficient optimisation of vehicle body treatments. These 

methodologies are essential to shorten development time and achieve high performance 
treatment configurations. They can be embedded into vehicle Computer Aided 
Engineering (CAE) design flow and can then be used in providing design and platform 
component sharing guidance information before prototype vehicles are available. 

15 

There are many methods known to the man skilled in the art to optimise damping 
treatments for vehicles. Nevertheless all these methods are based on standard Design 
of Experiment Techniques and are commonly gradient-based methods. None of them 
can efficiently and reliably handle a higher number of optimisation variables. 

20 

From EP 1 177 950 it is known how to optimise the vibration damping. This method uses 
a vibration damping material which is applied only in areas of maximum vibrational 
response after having determined maximum points of vibration response. The vibrational 
response of each section of the vehicle components is scanned after exciting some 

25 vehicle component from one point The disadvantage of this method is to be seen in the 
use of an experimental methodology for determining the vibration response. With the 
known method the damping material has to be applied on the areas with maximum 
vibrations, while the present invention looks to areas where the reduction of vibration 
shows a maximum overall effect, which is different than looking at the magnitude of the 

30 original vibrations. 

It is the aim of present invention to achieve a stable method for automatic optimisation ol 
vibration damping treatments for vehicles, which can efficiently and reliable handle a 
large number of variables. 

35 

This is achieved by a method comprising the features of claim 1 , and in particular by an 
optimisation tool named Genetic Optimisation for Lightweight Damping (GOLD). The 
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foundations of this tool are based upon Genetic Algorithms. The user can take into 
account many different damping material types and damping material thicknesses, 
distributed on metallic body panels of different thicknesses as well as different 
temperature areas the damping package. The damping package optimisation process \i 
5 automatically performed and controlled by the GOLD software. The user can 

interactively control the ongoing optimisation process and customise it by directly pre- 
setting special desired objectives and/or constraints, in terms of weight and NVH- 
performance (vibration or acoustic pressure), with simple interface commands. J 

10 Prsent invention, named GOLD, is a new too! for the automatic optimisation of vehicle 
damping treatments. It is based upon Genetic Algorithms, which are recognised as the f 
most powerful methodology to handle complex optimisation problems with very large 
numbers of possible discontinuous variables and exploits Finite Element (FE) and 
Boundary Element (BE) simulation techniques. The application field of this tool typical!: 

15 lies in the tow-middle frequency range. The possible variables taken into account for th^ 
damping treatment optimisation can be: 
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panel material; 
panel thickness; 
damping treatment type; 
damping treatment thickness; 

panel area temperature; for instance, higher temperatures can be assigned to 
body panels in the tunnel or dash areas or a temperature distribution map 
derived from a thermal camera can be taken into account; 
damping treatment local distribution. 
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In practice, GOLD can take into account a very high number of optimisation variables, 
which in theory would lead to a huge number of possible damping treatment 
configurations (typically billions) and efficiently handle optimisation problems practically] 
impossible to be solved with standard Design of Experiments (DOE) techniques and 
gradient-based methods. The necessary starting input data for the optimisation are the| 
following: 
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FE model of the structure (vehicle body) on which damping has to be optimised! 
In case an acoustic target is required, BE (or FE) model of the passenger 
compartment; 



- Material parameters of all the possible damping treatments used in the 
optimisation; 

- Definition of the damping patches, i.e. the possible areas on the vehicle body 
panels where damping can be applied; 

5 - rf required by the user, more additional constraints in terms of weight and NVH 
performance (vibration or acoustic pressure) can be defined. 

In the beginning of the optimisation process, the total panel area where the designer 
decides damping can be applied is subdivided into a set of sub-areas, named damping 

10 patches. A damping patch is just a possible pad of damping treatment. "Possible" mean? 
that the patch can either be treated or not GOLD can always take into account the 
possibility of leaving the patch bare. Each damping patch can cany a different 
combination of the above optimisation variables. Depending on the assembly or 
configuration of the damping patches, different damping packages are found. Damping 

is treatments are included in the vehicle FE model by means of the Emerald methodology, 
starting from measured frequency and temperature dependent material properties. 
However GOLD may be linked to any damping simulation methodology. 

Following a typically genetic evolution flow, the optimisation is performed in a cascade o 

20 selective 'iterations (generations). One damping package, named individual and coded 
as a binary string, is a possible solution of the optimisation problem. Each individual 
(binary string) corresponds to a treatment configuration of the body panels, i.e. a specific 
spatial distribution of selected damping treatments on the damping patches. At each 
generation step, the best performing individuals according to a defined objective functior 

25 (OF) are selected. The group of individuals is then able to generate new individuals by 
means of genetic statistical operators embedded into GOLD, like, for instance, 
crossover (exchange of some bit sequences) between the selected individuals or 
mutation (bit change for an individual). According to a statistical selection, the best 
performing individuals have the highest chance to reproduce, while the worst performing 

30 individuals have a low reproduction probability or may even be discarded and replaced 
with new ones. The initial population is randomly selected and its individuals are very 
different. This means global exploration of the optimisation domain takes place in order 
to achieve the global OF maximum; unlike gradient-based optimisation techniques, the 
optimisation performed by GOLD is not stopped when a local maximum is found. The 

35 vibration and acoustical response calculation of an individual is automatically run by 
GOLD in case a potential optimum individual is found. In the end of the optimisation 
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process, a group (population) of best individuals is kept in memory and the optimum 
damping package can be selected. 

The user has the possibility of monitoring the whole optimisation process on a visual 
5 user-friendly control panel where the state of the optimisation is shown in real time with 
its main parameters, iteration by iteration. The explored solution space (he- the domain 
containing all the possible solutions) typically has a dimension of m N , where 

- N is the number of possible damping patches 
10 - m is the number of possible treatment solutions 

The final optimisation target is typically the reduction of the following quantities: 

- the damping package weight 

. the vibration response with respect to frequency 
15 - the acoustic response (Sound Pressure Level: SPL) with respect to frequency 

The evaluation of the second quantity needs a structural FE run (i.e. Nastran®) 

The evaluation of the last quantity, whose procedure is described in the following lines. 
20 needs a structural FE run plus an acoustic BE or FE computation. 

In general, transfer functions from attachment points on the car body to the interior are 
usually given as targets for the body acoustic performance. The transfer functions are 
typically in terms of p/F y where p is the sound pressure at the passengers' ear locations 

25 and Fj is the force applied at the engine or suspension attachment points, in one 

direction. The p/F, transfer functions are calculated with a FE/BE uncoupled approach, 
which means that the fluid loading on the structure is neglected. The methodology can 
be divided into three main steps. In the first one, FE analysis is applied to calculate the 
vibration velocities v k at the nodes of the structure surrounding the cavity, due to the 

30 excitation forces F/. In the second phase, the BE analysis of the acoustic cavity is 
performed in i.e. Sysnoise® by placing a volumetric velocity unit source at the point 
where the p/F/will be calculated. The pressure at the nodes of the BE mesh is taken a 
output Thus, the plv k transfer functions are calculated reciprocally. In the third phase 
the velocities from the first step and the p/v* from the second step are combined to 

35 calculate the p/Fj amplitude and phase. 
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This technique provides some advantages in comparison to standard FE tools doing 
coupled vibro-acoustic calculations (like Akusmod®). A first big advantage is that the 
first two steps are independent If only the structure is modified, but the geometry of the 
passenger's compartment stays the same, the second step does not need to be 
repeated and vice versa. A second remarkable advantage is that the inner absorption of 
the panels in the passenger compartment can be taken into account locally, rather than 
through modal acoustic damping. Lastly, the contribution of the vibration body panels to 
the sound pressure can be calculated very easily. The main approximation is that no 
two-directional coupling is taken into account. In addition it can be that the BE element 
computation is time consuming. However, it should also be considered that the BE 
calculation is done only once and, if required, it can be replaced with an acoustic FE 
calculation. 

The optimisation process can be easily customised, with the possibility of choosing a set 
of pre-defined targets. The user can directly modify any of the above optimisation 
targets separately or choose any combination of them just by changing some command 
Ones contained into a simple text interface file. It is also possible to give priority to a 
special target instead of another. This enables the achievement of the best compromise 
between weight reduction and NVH performance improvement (either structural 
vibration or SPL), according to the particular constraints the user wants to apply. 

Two different algorithms are embedded into the GOLD optimisation toolbox: 

. The Standard Genetic Algorithm works on a relatively large population and 

converges with a more regular slope to the best solution. 
. The Micro Genetic Algorithm allows many generations on a much smaller 

population. It converges fast into a domain containing good solutions, then the 

slope of the OF curve decreases considerably. 

The efficiency and versatility of GOLD makes it easy to be applied to a wide range of 
damping optimisation problems involving simple test cases, automotive components or 
vehicle bodies. In the following a few typical examples are discussed: 

Damping optimisation on flat rectangular plate 

The very first time, the GOLD optimisation toolbox was applied to a relatively simple test 
case. The test structure was a free-free steel flat rectangular plate, size 800 by 500 mm 
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and 2 mm thick. A FE model of the plate was made first. Then the Mean Square Velocity 
(MSV) response of the plate carrying a given reference damping treatment configuration 
was calculated. 14 output points were randomly spread over the plate surface, after 
assessing that this was enough to well characterise the MSV response of the test 
5 structure. The excitation was a unit point force applied in a comer and perpendicular to 
the plate surface. The reference damping treatment was a 2 mm thick layer of damping 
material on the whole plate surface. After that, the surface of the plate was subdivided 
into 16 possible rectangular damping patches, all with the same surface size. 

10 GOLD was then used to perform a multi-objective optimisation on the plate and to find a 
new predicted optimised damping package giving lower MSV response levels and at 
least 20% lower weight than the reference damping treatment. Temperature was not 
considered as a variable in this case and was fixed at standard room temperature. For 
each patch, GOLD could choose among four different configurations: bare, 2 mm, 4 mm 

is and 6 mm damping. The genetic algorithm optimisation by GOLD could achieve an 
optimum distribution of damping patches with two different thicknesses on the plate: 2 
and 6 mm damping. GOLD managed to achieve significant improvements in the 
vibration behaviour of the plate with a contemporary 23% damping weight reduction in 
comparison with the initial treatment. The optimum damping package predicted by 

20 GOLD was then experimentally built up and measured. In the experimental set-up the 
plate was hung in order to reproduce free-free boundary conditions and was excited with 
white noise by a shaker in the same location and direction as the FE model, while small 
low-mass accelerometers were used for the output signal acquisition. 

25 The improvement in the MSV response predicted by the optimisation for the lighter 
optimised damping package could be experimentally verified, i.e. by comparing the 
values of the FE model and experimental set-up of the plate carrying the optimised 
damping package, together with the simulated and measured MSV responses. It can be 
seen that the experiments clearly confirmed that the prediction of MSV-reduction made 

30 by GOLD was correct 

Damping optimisation on vehicle component 

This GOLD application case on a vehicle component was the next step after the first 
35 validation case on the flat plate. The component was a part of a real vehicle floor with 
simplified boundary conditions. Its vibration response was firstly calculated taking into 
account the original damping treatment Comparative measurements of the untreated 



component for the same boundary and excitation conditions showed that the simulation 
was well able to reproduce the vibration patterns. In particular, it was important to 
assess that the areas with high vibration levels were the same in measurement and 
simulation, even though the levels did not perfectly match. GOLD was then used to 
5 perform the damping treatment optimisation on the component and find a new predicted 
optimised damping package giving improved vibration results, i.e. lower MSV levels, 
than the original reference damping treatment, keeping the same damping mass. The 
possible damping areas on the floor (areas which are possible candidates to apply 
damping treatment) were subdivided into 15 potential damping patches. For each 
10 allowed configuration of metal sheet thickness, damping type and thickness, a set of 
different input material properties for GOLD was created. Each patch could be left bare 
or carry three possible treatment configurations, which means that the total number of 
possible solutions of the optimisation problem was 4 1S »1.1 -10 9 . 

i s In order to reach convergence of the optimisation on reasonable calculation time, the 
Micro Genetic Algorithm was applied first on a reduced population of five individuals. 
After a sufficient number of generals, when an asymptotic behaviour of the Objective 
Function was reached, the Micro Algorithm was stopped, after about 1 hour. The 
population of individuals obtained at this step was used as input to the Standard Genetic 

20 Algorithm, which was run with a population of 30 individuals and 40 generations. After 
the GOLD run, which in total lasted just 10 hours, the distribution and mean square 
velocity response of the optimised solution could be examined. The optimum distribution 
consists of 11 treated patches out of the 15 possible candidates. Significant 
improvements in the mean square velocity vibration response could be predicted. After 

25 that, the predicted optimised damping package was manufactured on the real vehicle 
component and an experimental verification performed. Even though in this case some 
discrepancies between the simulation and experimental untreated component vibration 
levels was found, probably due to inaccurate modelling of the structure, the experiments 
confirm that the predicted damping package by GOLD performs globally better than the 

30 original reference package. 

The potential of GOLD for the spatial optimisation of damping treatments was also 
tested against a non-automated iteration of the "trial-and-error* kind. The "manuaF 
optimisation loops, mainly based on the experience of the user, allowed to find area and 
35 thickness distributions of the treatment with a little better mean square velocity response 
than the initial damping package. However, after a few days of calculations, the manual 
optimisation could not lead to any further vibration response improvement and still the 



optimised response was worse than GOLD. This can give a rough idea of the high 
efficiency of the automatic GOLD optimisation process compared to standard man- 
driven optimisation techniques. 

5 

Damping optimisation on vehicle body panels 

A real vehicle application example of GOLD is a damping optimisation performed on an 
upper class limousine car body. On this vehicle model a multi-objective damping 
10 treatment optimisation was performed in order to improve the vibration behaviour on the 
front floor while reducing the damping treatment mass by at least 30% with respect to 
the original damping package. 

Calculations on the whole car body and on the floor, using a substructuring technique to 
is model the remaining car body parts, were performed. In this way, more realistic 

boundary conditions can be applied to the panels (modal constraints at the nodes). In 
this particular test application a vertical unit displacement excitation was used. 

Firstly, the FE model of the panels without any damping treatment was validated again 
20 a set of measurements and the general good correspondence of the simulated panel 
vibration responses with the measured ones was positively verified. After this initial 
validation step, a FE calculation with the original damping treatment of the vehicle was 
performed and the MSV response extracted and used as the reference case which the 
optimisation had to improve. 

25 

The damping areas (areas which were possible candidates to apply damping treatment 
were divided into 17 potential damping patches. During the optimisation, each patch is 
allowed either to stay bare or to cany one out of seven possible different damping 
treatments (different materials and thicknesses). 

30 

The metal sheet thicknesses of the floor panels were given and different temperature 
areas for the dash, tunnel, front and rear floor regions were taken into consideration. In 
general, damping materials are strongly temperature dependent, thus the same 
treatment laying over panel areas with different temperatures has different mechanical 
35 properties: this behaviour is correctly taken into account by GOLD. In feet, it is possible 
to consider the same patch physically consisting of several sections having different 
thickness and/or temperature. 



A set of different input material properties for GOLD were created per each allowed 
configuration of: 

- metal sheet thickness 
s - damping thickness 

- damping material 
-temperature 

The optimisation domain given by the 8 possible treatment configurations with 17 
10 patches is made by 8 17 » 2.25 • 1 0 15 possible solutions of the optimisation problem. 

Of course, this figure includes all authorised configurations for the optimisation , some c 
which can be unrealistic, like for example "all 17 patches bare" or "ail patches treated 
with the highest treatment thickness". 



is In order to reach a good convergence of the optimisation in reasonable calculation time, 
the Micro Genetic Algorithm was applied first on a reduced population of five individuals 
with a high number of possible generations allowed. After the 67* generation a flattenini 
of the slope of the Objective Function versus the number of performed generations was 
remarked, so the Micro Algorithm was stopped at that generation, after about 30 

20 minutes. The population of individuals obtained with these steps was used as input to 
the Standard Genetic Algorithm, which was run with a population of 60 individuals and 
120 generations. 

After the GOLD run. the best solutions according to the optimisation criteria were listed 
25 and the optimised solution could be extracted. That solution consists of seven treated 
patches (out of the 17 possible candidates). Significant improvements of the MSV 
average and peak levels were achieved, together with 33% weight reduction of the floor 
damping treatment with respect to the original treatment When comparing the initial anc 
optimised MSV curves, three frequency domains with different efficiency of the damping 
30 treatment can be identified: 



below 100Hz, no local panel modes exist, so there is in practice no vibration 
reduction / damping effect of the treatment and thus there can be no improvement 
anyhow due to damping application 
35 - between 1 00Hz and 1 75Hz a moderate effect due to damping is visible 

above 175Hz, important improvements are achieved, especially near the mode 
peaks (up to 5 and more dB) 
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With the described strategy, a huge optimisation space of 8 17 possibilities could be 
efficiently explored performing 173 Nastran® runs: this meant in total 33 hours of 
elapsed calculation time on a standard Unix workstation. The gains in the number of 

s explored solutions and calculation time are enormous with respect to standard DOE 
techniques. In order to verify the prediction, together with the above GOLD optimisation, 
also an experimental optimisation of the damping package was performed by the 
Diamonds methodology on the real vehicle. The experimental optimisation process led 
to select the same damping material on the same vehicle panels with almost the same 

10 topological distribution of the treatment. However, a prototype of the vehicle is needed to 
carry Out the hybrid experimental-numerial Diamonds optimisation, while GOLD is a 
purely numerical procedure, which can be applied before any prototye has been made 
available. 

is The advantages of present invention are obvious for the man skilled in the art. GOLD is 
a new optimisation tool developed by applicant, based on a genetic algorithms and able 
to efficiently predict the optimum damping package on vehicle body panels in terms of 
materials, thickness and local damping distribution. GOLD allows the efficient 
exploration of very large solution domains and has the possibility of taking into account 

20 high numbers of variables, even in full vehicle computations (metal sheet and clamping 
treatment type, thickness, temperature and distribution). It has been successfully tested 
in different test cases, from simple plates to full vehicle applications. With its wide choice^ 
of algorithms and customisation functionality, it allows the definition of a wide range of 
optimisation strategies, according to the specific user needs (weight reduction, vibration 

25 reduction, improvement of the acoustics, improvement for a specific frequency range, 
etc.). The GOLD has an open architecture, which makes it easy to modify and link to 
any simulation methodology. It is compatible, for instance, with all MSC Nastran® 
features, in particular sub-structuring techniques (superelements) and modal 
superposition calculations. 

30 

GOLD helps the achievement of the optimum in the design of vehicle damping packaged 
by means of FE simulation. It is based on generic algorithms, which are recognised as 
the most powerful methodology to handle complex optimisation problems with a very 
large number of variables (possibly discontinuous). The optimisation is performed in a 
35 cascade of successive iterations. In each iteration a set of solutions is explored and the 
best are selected on the basis of certain given constraints: lower vibration, lighter 
weight, tower SPL. Each time a new iteration is performed, the population evolves for j 

sinn4s.<i 



fa 



the better the best solutions are kept in the selection process and can generate new 
solutions, while the worst are discarded and replaced with the new ones. The whole 
iteration and generation process is automatically performed and controlled by special 
operators embedded in the GOLD software, which is written in MatLab® software 
s package. 

The user is free to customise the damping package optimisation by directly pre-setting 
the optimisation constraints with simple interface commands. The user also has the 
possibility of monitoring the optimisation process on a visual control panel, where the 
10 state of the optimisation flow is shown in real time, iteration by iteration. One of the main 
features of GOLD is that the number of input variables can be very high. The user can 
take into account many different materials, distributions and thicknesses of metal and 
damping, as well as different temperature areas in the damping package. 
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1 . Optimising method for vibration damping treatments of vehicle body panels 
characterised by the following input steps to create input variables: 

generating a structural finite element (FE) model of the vehicle body 
respectively the panel(s) (dash, floor, tunnel) on which damping has to be 
optimised and/or generating a boundary element (BE) model of the passenger 
compartment (in case that an acoustic target/SPL is required); 

defining damping patches potentially subject to possible treatment; 



defining a plurality (N) of possible damping treatments (including no treatment) 
and determining its material parameters / characteristic data as well as the 
is body panel material parameters / characteristic data; 

if required, defining more additional constraints or objectives in terms of weight 
and noise vibration harsh (NVH) performance (vibration and acoustic pressure) 

20 - and further comprising the following computing steps: 



computing the equivalent material damping parameters / properties (m) (i.e. 
overall thickness, weight, porosity, bending-stHfhess, elongation stiffness, 
bending loss factor, elongation loss factor, visco-elasticrty, temperature, etc.) bj 
25 multi-layer simulation (for instance by EMERALD 9 ) from these material 

parameters for each combination of any of the plurality of the possible damping 
treatments with any of the vehicle body panel parameters, including the 
temperature and frequency dependency of all materials involved; 



30 - running a FEM simulation (i.e. a finite element model simulation, such as 
NASTRAN®) for a reference configuration in order to calculate the dynamic 
response (NVH response transfer function), in particular the vibration behavioui 
of the vehicle on which damping has to be optimised with respect to exitation, 
architecture/structure and used materials; 

35 

applying a genetic algorithm to the above input variables, which algorithm 
performs a selective iteration by: 
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a) generating a pool of individuals (damping packages/treatment configuration, 
coded by a binary string) from the input variables (equivalent material damping 
properties, spatial distribution, thickness, weight, etc. of the binary string); 

b) mutating (bit change for an individual)and/or crossing (exchange of bit 
sequences) randomly selected individuals/genes of this pool by means of 
genetic statistic operators to generate a new generation of individuals/genes; 

c) selecting each individual of the new generation by means of a statistical 
selection according to a defined objective function (OF), i.e. calculating the 
value of a predetermined fitness/priority criterium/predefined targets (lower 
weight, lower vibration, lower sound pressure level, lower cost, etc/objective 
function); 

d) correlating individual's chance of mutating and/or crossing with their 
performance with respect to the objective function; 

e) mutating and/or crossing these chance-correlated individuals by means of 
genetic statistic operators to generate a next new generation of individuals; 

f) iterating steps c), d) and e) until a predetermined flattening of the slope of the 
OF versus the number of performed generations is achieved, which leads to a 
set of optimised vibration damping configurations. 
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